Nematicity driven by hybridization in iron-based superconductors 
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In this paper we study an effective model for the normal state of iron-based superconductors. It 
has separate, but interacting itinerant and localized degrees of freedom, originating from the dxz 
and dyz, and from dxy iron orbitals respectively. At low temperatures, below a mean-field phase 
transition, these different states condense together in an excitonic order parameter. We show that at 
even lower temperature, after another phase transition, this ordered state can spontaneously break 
the C4 lattice symmetry and become nematic. We propose this mechanism as an explanation of the 
tendency towards nematicity observed in several iron-based compounds. 



Introduction. The discovery of high-temperature su- 
perconductivity in iron-based material^iM^ is one of the 
most exciting recent developments in physics. It is almost 
certain that the origin of superconductivity is unconven- 
tional (i.e. driven by electron-electron interactions), and 
is very likely that the order parameter is unconventional 
as weli'' -^L The normal state properties of these mate- 
rials, in contrast, at first seemed rather unremarkable ~ 
a bad metal behavior at high temperatures, which can 
be followed by structural transition and antiferromag- 
netic phase at low temperatures. This simple picture 
was considerably complicated by the growing evidence of 
nematic fluctuationJ^Hs] or even long-range ordei'SJ present 
in some compounds at temperatures well above the an- 
tiferromagnetic and structural transitions. Furthermore, 
this nematicity seems to be of purely electronic (rather 
than structural) origirP^l, Understanding this state and 
its connection to the other phases is an important mile- 
stone in the study of these materials. 

Several scenarios have been proposed to explain this 
tendency towards ncmaticity^^. Some suggest that the 
nematicity is generated by the spin fluctuations, and is 
a precursor of the incipient antiferromagnetic stat^^^HSl 
Another class of models emphasises the propensity to- 
wards orbital order of the iron Sd states in different 
environment^i^lliSl There is also a model which considers 
the nematic state a consequence of a more complicated 
"hidden" ordeM 

In this paper we suggest an alternative scenario for 
the nematic instability. We introduce a simple effective 
model of the normal state of these materials, that as- 
sumes the simultaneous existence of localized and itin- 
erant electronic states (originating in the iron d^y and 
dxz/dyz orbital states respectively). Such models have 
already been considered in the context of iron-based 
superconductor^^HHUl concentrate, however, on the 
mixing of these two types of states, driven by hybridiza- 
tion and interactions. We show that interactions create 
a tendency for an excitonic instability of the localised 
state with more itinerant bands on the neighbors: at 
the mean field level this generates a hybridization gap 
with the symmetric dxz/dyz combination which physi- 
cally represents a crossover to a screened orbital at low 



energies with no symmetry change. At lower temper- 
atures, we find that such order parameter can sponta- 
neously break the C4 lattice symmetry, and thus de- 
scribe a nematic state. This scenario emphasizes the 
coupling of dxy with dxz and dyz states, and naturally 
explains the recent an gle-res olved photoemission spec- 
troscopy (ARPES) datcP^l2SI, in addition, a hybridiza- 
tion gap opens at the transitio n, wh ich may explain some 
of the ARPE^, conduct ivitji^IlSSI^ opticaP^ and point- 
contact spectroscopji^Sl data. 

Effective model. To arrive at the effective model we 
start from a realistic tight-binding Hamiltonian for the 
three t2g iron d orbitals (the Cg orbitals do not contribute 
significantly to the density of states close to the Fermi 
level and we neglect them), which can be written aPl' 

V, = Tio + Hint- (1) 

We use the one-iron-atom unit cell axes as our coordi- 
nate system (and note that dxz and dyz orbitals are de- 
generate). Ho contains terms like J2ij,ai3atai3dlacrdj,i3a, 
with tunnelling amplitudes (where i, j and a, P are 
site and orbital indices, respectively) which include the 
effects of the pnictide/chalcogenide atoms. We assume 
that all orbitals are close to half-filling. T-Lmt represents 
the SU (2) invariant on-site interactions, which includes 
the intra- and inter-orbital density-density interaction 
and the Hund's couplin^^. 

One key insight from the dynamical mean-field the- 
ory (DMFT) calculations is that in some of these com- 
pounds the dxy orbitals experience the so-called "kinetic 
frustration" due to the destructive interference of several 
tunneling paths for the dxy electron^2ll. This suppresses 
the di^xy ^ dj^xy hopping, and makes tl^y^xy significantly 
smaller than the other tunneling amplitudes. Thus the 
dxy electrons are much more localized then the rest, and 
in some cases may even undergo an orbitally-selective 
Mott transitioiP^. 

Based on this we suggest a simplified effective Hamil- 
tonian, which nevertheless preserves some of the key fea- 
tures of the realistic model described by Eq. ([!]). It con- 
sists of four parts; 

^e// = T~{-c + 'Hd + 'Hrnix + l~ires- (2) 
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The first term describes bands of itinerant electrons (de- 
noted by c), which originate from the dxz/dyz sector of 
the model: 

V-c^Yl ek.p4,A.<TCk,M^, (3) 

where /i is a band index. The c states are product of the 
mixing of the dxz and dyz orbitals. Even after including 
some interaction effects, we assume that the renormalized 
band structure is similar to that of Ref. 1351 Since we 
are interested in the low-energy physics we treat the four 
surfaces crossing the Fermi level as separate bands. Thus 
(X runs from 1 to 4, with 1 and 2 denoting the hole pockets 
around T and (tt, tt) points, and 3 and 4 the two electron 
pockets around (0,7r) and (7r,0) points (in the so-called 
unfolded Brillouin zone). 

The second term in "He/ / contains localized states (de- 
noted by d) on a periodic lattice, with large on-site re- 
pulsion: 

■Hd = ^ Eodl^di^ + Ui^dani^da' , (4) 

where rii^a- is the density of d electrons with spin a on 
site i. It describes the narrow, strongly correlated dxy 
band, present in some of the pnictides and chalcogenides. 
The dxy states are taken as completely localized, and we 
assume Eq to be below the Fermi level of the itinerant 
bands. 

Hmix describes the mixing between the two types of 
fermions: 

'H,nix = V^{dl„Ci^^„ + cl^^di„) + W^ni^dni,c,{5) 

% i 

where c| creates an electron in a Wannier state of band 
/i. One part of it is hybridization, originating from the 
xy O xz and xy yz hopping terms in T-Lq^ which are 
not necessarily small. Note that we are considering only 
real on-site V (since the c states are a superposition of dxz 
and dyz orbitals from different sites the on-site hybridiza- 
tion is allowed). More realistic V has to be complex and 
have a long-range part. The other - four-fermionic ~ 
terms come from the interorbital interactions, and we as- 
sumed a simple density- density form for them. All terms 
in "Hint which cannot be cast in such form we hide in 
T-Lres ("res" for residual). We neglect these interactions 
for the moment, but they can play an important role. 

As a result of this simplifications we arrive at a model 
of periodic lattice of localized states embedded in a four- 
band "sea" of itinerant electrons. Similar models have 
already been discussed in relation with the magnetism 
of the iron-based materiald^2tt23j^ however, the term T-Lmix 
has usually been omitted. As we will show, it has impor- 
tant consequences and can lead to interesting effects in 
the paramagnetic state of these materials. 

Mean-field approximation. We use the Hamiltonian of 
Eq. ([2]) to study the normal state of iron pnictides and 
chalcogenides, well above the superconducting and the 



magnetic critical temperatures. Unfortunately, even this 
reduced T-Lef / cannot be solved exactly. 

We consider the limit of very large C/, in which the 
doubly-occupied d states are effectively forbidden. To 
deal with this constraint we employ the slave-bosons 
techniqu^^ by the transformation d^^ — > dicb\. There 
is also an on-site constraint ^„d\^dicr -\- b\bi = 1, 
which can be enforced by a Lagrange-multiplier term 
Ai(rtj + b\bi — 1) in the Hamiltonian. 

Without the interaction term in Umix the effective 
Hamiltonian in this limit is very similar to the Periodic 
Anderson Model (PAM), studied extensively in relation 
to heavy- fermion materials'^^. The mean- field treatment 
of the problem leads to qualitatively correct result for 
these compounds (with some important caveat^^sESI) 
Because of this analogy we expect at low temperatures 
the d electrons in our model to delocalize due to the 
mixing with the itinerant bands through V (there are 
evidence of such crossover in the normal state of iron 
chalcogenidePi'^ . This delocalization is described on a 
mean-field level by a phase transition, and the appear- 
ance a non-vanishing expectation value of (b) = b, and 
can significantly change the electronic spectrum by open- 
ing a hybridization gap. 

To demonstrate this we introduce two mean fields - 
b and A. b describes coherent tunnelling of d electrons 
between neighbouring sites and A enforces the no-double 
occupancy constrain (but only on average). However, we 
still have to deal with the four-fermion terms in Umix- 
We can decouple them by introducing additional mean 
fields in the particle-hole channel (j)if^ = (dl^Cia). As- 
suming spatially uniform case we use the replacement 

W J2 ^d^^^c,r -W Cjy^cl^j,, - . (6) 

It is important to note that (j)^ should be thought of 
as components of a four-dimensional vector, rather than 
four independent order parameters. The reason is the 
obvious constraint that one d state can only condense to- 
gether with a single c state or a particular linear combi- 
nation of the different c^i (which one is determined by the 
particulars of the model). However, the remaining three 
orthogonal linear combinations cannot condense for the 
lack of available d electrons (similar physics is behind the 
so-called "hidden metalicity" in the three-band model of 
Ref. HSI). 

Hybridization gap. We first consider a fully symmet- 
ric version of our model, with two identical circular hole 
bands around F and (tt, tt) points, and two identical cir- 
cular electron bands around (tt, 0) and (0, tt) points. For 
simplicity we assume that 4>i = 4>2 = 0, i.e. no hybridiza- 
tion of the d electrons with the hole bands (adding these 
fields is straightforward). The fact that the two electron 
bands are identical seems to imply that the order param- 
eter can be freely rotated in the (3, 4) subspace. However, 
as shown in the next section, we have to take symmet- 
ric combination of (f>i and 02, or using vector notation 
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(0, 0, (/), (/)). Note that this order parameter preserves the 
underlying lattice C4 symmetry (due to the 3 O 4 i.e. 
X y symmetry of the order parameter). Also, as al- 
ready explained, it leaves the other three "bands" - ci, 
C2 and the antisymmetric combination of C3 and C4 - de- 
coupled on mean-field level. 

Thus, with the help of the static and uniform fields 
A, b and <j) we have reduced the effective problem to a 
solvable single-particle model 

Umf = ^ ekCk,o-Ck,o- + '^{Eo + \)d\^d.,„ - A + 

kcr i(7 

^(A62 + W<t>^) + J2 Vbdlc,^, - J2 W<Pdlc., (7) 

i ia ia 

(plus three free-fermion Hamiltonians), where we have 
dropped the tildes of the slave-fermion operators. The 
self-consistency equations for A, h and (j) are 



VT Gdc + \h = Q, 

k.fi.ujn.a 



= 0, 



(8) 
(9) 
(10) 



where we have defined the following Green's functions: 



Gdc — 
Gd^ 



V 



[iOJn - ^k)[ii^n - e.d) - ' 
(iw„ - efe)(iw„ ~ td)~y'^' 
(ia;„ - ek)iiuJn - ed) - 



(11) 
(12) 

(13) 



with ed — Eo + X and V — Vb ~ Wcj). We solve these 
equations to obtain the behavior of the mean fields with 
temperature, and show the result on Fig. [l] As can 
be seen, there is a phase transition at temperature Tc- 
Above it we have 6 = 0, = and A = 1 - the itinerant 
and the localized states are decoupled. Below Tc b and (p 
simultaneously develop non-zero expectation value, and 
A becomes bigger than 1; this describes the delocaliza- 
tion of the d^y states through coupling with the dxz/dyz 
bands, and the opening of a hybridization gap V in the 
electronic spectrum (also, note that (j) is negative - it 
increases the gap, and even dominates it at low temper- 
atures). This phase transition, however, is an artefact of 
the mean field treatment ~ finite b violates the local gauge 
symmetry of the model. Formally, infrared fluctuations 
destroy the order and restore the symmetrjl^SttiO] Never- 
theless, as in the case of heavy-fermions we can take the 
mean field transition as an indication of a crossover in 
the real system. 

Breaking of the C4 symmetry. To see how the sym- 
metry between and 4)4 appears, and how it can be 
broken we consider a Ginzburg-Landau (GL) theory for 




FIG. 1. Plot of the solutions for (purple), b (green) and A 
(red dashed) of the self-consistency equations'*^- that show the 
mean- field phase transition at T^. The delocalization of the d 
electrons is described by the appearance of finite negative </!> 
and positive b. The inset shows a zoom of the Tc region. 



the fields b and 0^ (again assuming = 02 = for 
simplicity). For the symmetric case of identical circular 
electron bands the GL free energy takes the form: 

J-Q = abb'' + I3bb^ + (a^|0|2 -f /3^|0|4) 

-H^7i6|0^|cos^^, (14) 

where is the phase of 0^. The coefficients in this 
(bosonic) theory can be obtained by integrating out the 
fermions, and depend implicitly on A. 

Now it is clear why we had to take the the symmetric 
combination of 03 and 04-6 acts as an external field 
in the (1, 1) direction of the (03,04) space. Thus the C4 
symmetry is preserved both above and below the tran- 
sition. The linear coupling also guarantees that b and 
become non-zero simultaneously, and justifies our choice 
of mean-field decoupling of the W term. 

In realistic models the symmetries of the band struc- 
ture are significantly reduced. First, the hole pockets 
have different kp (thus 0i and 02, if present, are not de- 
generate). Second, the electron pockets become ellipses 
(but are still identical after translation and tt/2 rotation). 
In terms of orbital content the Ci and C2 states are still 
equal mixtures of dxz and dyz (averaged over their re- 
spective Fermi surfaces). In contrast, the C3 and C4 are 
dominated by either one or the other orbital. Impor- 
tantly, this means that any order parameter that is 
not 3 -H> 4 symmetric will simultaneously break the C4 
symmetry of the lattice, and induce orbital order. 

It is easy to see that the ellipticity of the bands leads 
to additional quartic term in the effective action: 



-F' = E^2i0.n0Mi 



(15) 



with ^ ^ v, through the standard four-fermion diagram 
on the left of Fig. [2] The value and the sign of the co- 
efficient depend on the ellipticity, but it can be shown 
by direct calculation that for parabolic bands with small 
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ellipticity 72 is positive. The physical origin of this cor- 
rection is the particle-hole asymmetry. 




FIG. 2. On the left is the Feynman diagram that contributes 
to 72 because of the ellipticity of the C3 and C4 bands. On 
the right is diagram of a process generated by the interband 
interactions. 

There is yet another source of quartic terms - the 
renormalized interactions in the d^z/dyz subspace hidden 
in 'Hres- In particular, c|^c^cJ,Cy term generate processes 
like the one shown on the right of Fig. [2] They also 
contribute to 72, and again it can be easily shown that 
their contribution is positive. Note that these are inter- 
band interactions terms, which prefer the delocalized d 
electrons to mix with C3 or C4; the (presumably larger) 
intraband terms participate in setting the overall scale 
the order parameter, but do not contribute to 72. 

If 72 is positive it penalizes the coexistence of 03 and 
(/>4, and favors an order parameter with only one non- 
zero component. The effect of this term, however, can 
overcome the linear coupling of ^3 and 04 to b only at 
a finite temperature below Tc; this is demonstrated on 
Fig. |3] in which we show solutions of the GL equations 
below Tc- We see that at temperatures sufficiently below 
Tc an asymmetric combination of 03 and 04 becomes a 
minimum of the free energy. This state breaks the C4 
symmetry of the underlying lattice and reduces it to C2. 

Thus in our model the nematic state emerges from a 
two-stage process - at the first stage the dxy electrons, 
localized at high temperatures, start to mix and form 
coherent bands with the itinerant electrons. The quasi- 
particles that emerge at this stage are strongly renormal- 
ized by local interactions like Hubbard's U and Hund's 
Jh- Note the parallels not only to the PAM, but also 
to the Zaanen-Sawatzky- Allen worlP^, in which c?-orbital 
electrons delocalize through coupling with valence bands. 
At a lower temperature, possibly far below the mean-field 
transition, the order parameter breaks the lattice rota- 
tion symmetry. This is the only true phase transition 
of the model, and it is driven by band-structure details 
and the interband interactions of the itinerant electrons. 
Those are relatively weak effects, compared to the lo- 
cal physics dominating at higher temperatures, and they 
can play a nontrivial role only because the first transition 



(actually crossover) has significantly reduced the effect of 
the strong local interactions on the low-energy physics. 

In this scenario the nematic state is not induced by 
proximity to the antiferromagnetic state, but is an inde- 
pendent instability. The presence of nematicity, however, 
can, by itself, enhance the magnetic instabilitjS3. Also, 
note the similarities of our model with the physics of 
some nickel-based compounds^^ 
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FIG. 3. Plot of b (green), 03 (red) and 04 (blue) that minimize 
the GL functional Below T « O.MTc the effect of positive 
72 overcomes the linear coupling with b and additional solu- 
tion appears, with 03 7^ 04. In the inset we show that the 
free energy of this solution (red) is below the free energy of 
the symmetric solution (purple). The new order parameter 
breaks the C4 symmetry. 

Conclusion. In this paper we introduced an effective 
model for the normal state of iron-based superconduc- 
tors. It has both itinerant and localized degrees of free- 
dom - the former originate from the d^z/dyz, and the 
latter from d^y iron orbitals. We studied this model on a 
mean field level and showed that at low temperatures 
the dxy states can effectively delocalize and condense 
together with the itinerant states in an excitonic order 
parameter. Because of the multiband character of the 
itinerant Fermi surface the coupling between the local- 
ized and itinerant electrons can naturally lead to another 
phase transition, and a nematic excitonic state. We pro- 
pose this mechanism as an explanation of the tendency 
towards nematicity observed in several iron-based com- 
pounds. 
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